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Abstract 

We consider the exact solution of a many-body problem of spin-s particles interacting 
through an arbitrary U(l) invariant factorizable 5"-matrix. The solution is based on a 
unified formulation of the quantum inverse scattering method for an arbitrary (2s + 1)- 
dimensional monodromy matrix. The respective eigenstates are shown to be given in 
terms of 2s creation fields by a general new recurrence relation. This allows us to derive 
the spectrum and the respective Bethe ansatz equations. 
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The existence of exactly solved models have been playing a fundamental role in our under- 
standing of interacting one-dimensional many-particle systems in distinct areas of theoretical 
physics m El El II]- An important class of models are those constituted by n s species of particles 
of identical mass that interact only through a sequence of two-body scattering processes. The 
two-body collision amplitudes S(pi, A^)^ 1 '^ 2 , where aj,bj = 1, . . .,n s represent internal quan- 
tum numbers, are responsible for the redistribution of momenta among the different types of 
particles [5]. In general, these theories lead us to deal with non-diagonal scattering matrices 
defined by, 

£12(^1,^2)= Si^^lltl^ax®^^ (1) 

ai ,02 ,61 ,62 

where e a ,b denotes n s x n s matrices having only non-null element with value 1 at row a and 
column b. 

The rapidity fij parameterizes the one-particle momenta p(fij) and energy e([ij) associated 
to the jth particle. We stress here that we will consider the most general situation in which 
the S-matrix (jJJ depends on both values of two independent p,\ and /12 scattering parameters. 
The scattering from the initial state (01,02) to the final state (61,62) may exchange quantum 
numbers and momenta but do not alter the values of the latter parameters. This type of 
factorizable theory has been first emerged in the study of the eigenspectrum problem associated 
to particles interacting through delta-function potential [HI IE] • It was then discovered that in 
order to assure exact solubility the S'-matrix (JTJ) needed to satisfy a necessary condition called 
Yang-Baxter equation |H] , 

where sum of repeated indices is assumed. Here we shall consider solutions of (J2J) that are 
almost unitary, namely 

for some arbitrary function p^fix,^). 
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Assuming the existence of such integrable systems, the next step is to quantize them to 
determine the respective eigenspectrum. As first argued by Sutherland [9j this can be ac- 
complished within the framework of an asymptotic Bethe ansatz approach. The basic idea 
is that in the regions of the configuration space where the distance between the particles is 
much larger than the interaction range, the wave-function can be considered as that of a su- 
perposition of plane waves. It can be constructed similarly as for the problem of particles with 
delta- interaction [Sj, by replacing the scattering matrix of the delta-potential by a general fac- 
torizable S'-matrix If we consider periodic boundary conditions on a ring of size L, the 
asymptotic one-particle momenta p(fJ.j) is required to satisfy the following eigenvalue equation, 

where N is the number of particles in the system and ~~\t denotes the set {/ii, /z 2 , • ■ • , Hn} 
of corresponding rapidities. The energy E of the system is given by adding the one-particle 
energy expressions, namely 

N 

£ = £e(^). (5) 

The function A(A, ~ft) represents the eigenvalues of an auxiliary [n s ] x [n s ] N operator 
T(A, ~]f) usually called transfer matrix. This auxiliary eigenvalue problem can be defined by, 

T(A, Jt) \j>) = Tr^ [T A (X, Jt)} \4>) = A(A, Jt) \j>) , (6) 

where the trace is taken over an auxiliary n s -dimensional space A = C ns . Furthermore, the 
monodromy operator T^(A, ~~jt) is related to the S'-matrix elements by the following ordered 
product P1IJEEI], 

Xa(K A* ) — San(K ^n)S_an-i(^, A*at-i) • • • Sai{K A*i)- ij) 

In order to make further progress it becomes crucial the exact solution of the eigenvalue 
problem ()6I7|) for arbitrary n s and two-body scattering amplitudes. This is indeed a tantalizing 
open problem, especially when a particular form of S(/ii, ^2)0 & is not specified. Indeed, most 
of the results concentrate on specific S-matrices such as those related to the the six-vertex 



model [12], to its higher spin descendents [TSj and those based on higher rank Lie algebras 
[LH Ho"! \W\ IT7| ITS] . In the latter case, some of the findings ^7J EI] are still in the form of 
conjectures for the transfer matrix eigenvalues and also the number n s actually encodes a 
variety of distinct conserved quantum numbers such as spin, color, flavor, etc. Consequently, 
the number of null scattering coefficients grows rapidly with n s due to the many possible 
underlying U(l) symmetries. 

In this work, however, we shall establish the essential tools to solve the eigenvalue problem 
()6l7j) when only a unique U(l) symmetry is present for arbitrary n s . This is the minimal 
continuous invariance one could request and our results will be valid for arbitrary factorizable 
S'-matrices satisfying such symmetry condition and the unitarity relation (J3J). More precisely, 
we are considering integrable models whose S'-matrices fulfill the property, 



5*12 (/^l, 1^2)1 S{ + 5*2 



0. 



where is the azimuthal component of spin-s operator associated to the jth particle such that 
s = (n s — 1)/2. Note that relation (JBJ means S(X, = unless the ice rule a± + a 2 = 61 + 62 

is satisfied which leads us to a total number of n s (2n s + l)/3 non-null amplitudes. 

The direct connection between the number of species and the spin of the particles makes 
these integrable models physically relevant. We shall tackle this problem by means of the 
quantum inverse scattering method fUJE]- We remark that, in general, there is no known 
recipe to perform this task, therefore our results can be considered as new developments in 
this framework. 

The most important quantity in this method are the monodromy matrix elements on the 
auxiliary space which here will be shortly denoted by T(A) M , p,q = l,...,n s . These are 
operators on the quantum space TI^Li ®C" S such that the diagonal entries define the transfer 
matrix eigenvalue problem (jUJ) while the off-diagonal ones will play the role of creation and 
annihilation fields. The commutation relations between such matrix elements are given with 
the help of the corresponding Yang-Baxter algebra, 

S(X, ^TaA^M = %MT a)d (X)S(X, vfl*. (9) 



In order that the model be soluble by means of an algebraic Bethe ansatz, it is fundamental 
the existence a vacuum state |0) such that the monodromy operator (JJJ) acts on it as a triangular 
matrix on the auxiliary space for arbitrary A. Thanks to the underlying U(l) symmetry, it is 
possible to build up this state by the tensor product of local vectors, 



TV 



io> = n® !*>,■> 



(10) 



3=1 



where |s) ■ denotes the jth spin-s highest state vector. 

It turns out that this is one possible eigenvector on which the matrix elements of T Ptq (\) 
operate as follows 



%M 10) 



N 



3=1 

for p < q 



\pq) for p < q 

where \pq) are non-null vectors representing excitations over the vacuum |0). 

From the above expressions we see that the fields T p g (A) for p > q acts as creation operators 
with respect to the reference state |0). Of particular importance are the operators 7j g (A) 
(q > 2) which satisfy the following property, 



N 



2i,*(A),5>i 

3=1 



(g-l)T llfl (A). 



(12) 



A direct consequence of (I12j) is that the fields 7i i9 (A) can be interpreted as raising operators 
associated to excitations over the ferromagnetic vacuum |0) with spin component s — q+1. It is 
therefore plausible to suppose that the eigenvectors = $ m (Ai, . . . , A m ) |0) of © should be 
constructed similarly as a Fock space having 2s possible distinct creation fields. Considering 
that both the total number of particles and spin are conserved quantities one expects that the 
eigenvectors structure should be as follows. The field (Ai) represents the one-particle state, 
the linear combination 7i ) 2(Ai)7i )2 (A2)+V'i(Ai, A2)7l 3^1)7^1^2), for some function A 2 ), 
the two-particle states and so forth. It turns out that the form of such linear combinations 



can be inferred on the basis of the commutation rules between the fields 7i i9 (A) derived from 
the quadratic algebra (JHj). In addition to that, we emphasize that much of the simplifications 
needed to construct suitable eigenstates are carried out only on basis of an extensive use of the 
Yang-Baxter (J2J) and the unitarity Q constraints between the scattering amplitudes. Omitting 
here the technicalities of these computations ^H] we find that the multi-particle states obey a 
2s-order recursion relation given by, 

min{2s,m} 

$m(Ai, ■ ■ ■ , A m ) = ijik (Ai, Xj 2 , . . . , \j k )T 1 i +k (\ 1 )§ m _ k (\ 2 , . . . , Aj 2 _i, Aj 2+ i, 

k=l 2<j 2 <---<j fc <m 2 

k k 

. . . , Xj 3 -i, \j a +i, ■ A Jfc _i, A ifc+ i, . . . , A m ) x J[ Ti^(\ jd ) J[ 

d=2 1=2 

x ft §H^S ft e(A*,A*), (13) 
where we identify $o = 1 an d function 6(A,/i) is, 

>3.1n/i \2.2 o/\ \2,2 /\ \3,1 



5(A, ^) 3 ;i ^(A,/x) 2 ; 2 - s(x,fi) 3 ' t l s(x, /i) 2 ; 2 

'S'(A,^)i;i5'(A,/i)3;i 



( A ^) = ! ~ aTT^T -■ ( 14 ) 



In the course of our analysis we also have made the natural hypothesis that the fields states 
()13j) are symmetric functions in all {Xj} variables. This exchange property between the variable 
Ai and A 2 can be used to determine if)k(\i, Xj 2 , . . . , Xj k ) in terms of the scattering amplitudes 
recursively. The simplest case being s = 1/2 where by fixing V'i(A) = 1 we find that ipi(X, fj) 
is. 



^i(A,/i) = -— — -3V. (15) 



Similar task for higher spin is in general more involved but it can be done for any specific 
value of the spin [T§] . As an extra example we present below the explicit expression for s = 3/2, 



?MA,/i,t) 



6(A, /x)^i(A, t)S(X, n){jS(X, 

s(x,fi)y 2 s(x,fi)lj- s(x^)lis(x,fi)tj 
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For s = 1/2 our result (jl3j) recovers the known algebraic Bethe states associated to the 
six- vertex model fOJ IH| while for s = 1 we have an extension of the Bethe states for nineteen- 
vertex models |2()j . In fact, apart from factorizability and unitarity, no other assumptions 
on the /S-matrices elements entering the eigenstates (j!3|) have been made. In addition, the 
generality of our recursive manner to generate the eigenstates for arbitrary spin-s is, as far 
as we know, a new progress in the quantum inverse scattering approach . It follows from 
commutation relation (fT2"J) and our construction (fT3|) the property, 



which corroborates the physical interpretation of $ m (Ai, . . . , A m ) |0) being multi-particle states 
over |0). 

We now turn our attention to determination of the eigenvalues A (A, ~ft). From © one has 
to carry on the diagonal fields T PtP (X) over the creation fields Ti )q (fi) (q > 2) that built the 
multi-particle states (jlHj) . This task is made by recasting the Yang-Baxter algebra (jHI) in the 
form of commutation rules between these fields. In general, convenient commutation rules do 
not follow immediately from Q and a two-step procedure is needed. To give an example of 
our approach let us denote by [/; k] the Ith row and the kth column of (jHJ). The appropriate 
commutation relations between the fields 2 (A) and T PtP (X) for 1 < p < n s are obtained by 
using the combination [I] (I — 1) *n s + 2]S l l +l'l(X, /i) — [I; Hris + ^S 1 ^ 1 (A, fi) of the entries of 0. 
The basic idea is to keep the diagonal operator T P ^(X) always in the right-hand side position in 
the commutation rules and such procedure can be implemented for all creation fields 7i i( j(//). 
The eigenvalues are easily collected by keeping only the first terms of the commutation rules 
among 7i g (/i) and 7^ jP (A) and after some cumbersome simplifications we find that the final 
result for A(A, Jt) is, 



A' 



£ Sj$ m (\i, . . . , A m ) |0> = (sN - m)$ m (A 1 , . . . , A m ) |0> , 



(17) 



A(X,jt) 



i=l 1=1 



A' 




]l k=2 i=l ' 1=1 



n s -l iV m 



+n^/*)isn 



i=l 1=1 



N m 



S(X,\{) 
S(X,Xi) 




(18) 
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where 

^k{X,V) - — — ' M , — Tk+T^ 

The remaining terms that are not proportional to the eigenvector can be canceled 
out by imposing further restriction on the rapidities {Xj}. These are auxiliary Bethe ansatz 

equations and in our case only a unique relation is needed to eliminated all such unwanted 
terms. It is given by, 

n V S(Xj, fii)l'l _ tt rv\ \ \^C\/> ^( Ai ' ^3)2,1 fon \ 

— — _ ^ ^{Aj, Aij— — r~2j , ^1,1 • l^UJ 

1=1 S\ A ji Hl)2,l i=l ^l A j5 A ij2,l '- , l A i5 A iJl,l 



Before proceeding we remark that our expressions (|18II20|) reproduce the eigenvalues and 
Bethe ansatz equations of a particular integrable model with higher spin obtained through 
fusion procedure of six- vertex S'-matrices [T5] . Even in this special case, one hopes that our 
general expression for the eigenstates (JT3J) could still be of utility as far as correlation functions 
are concerned [2*Tj . 

At this point we have been able to derive that basic ingredients to obtain the quantization 
rule for the one-particle momenta p{fJ>j). This follows directly from ()4I18|) and the unitarity 
property of the S-matrix (JHJ), i.e. the condition S(X, X) c ^ d b ~ 5 a ,dSb, c - The final result is, 

e ^ = n 5^,^)55 fillr^ (21) 

i=l 1=1 b\ A h Hj)2,l 

The results (|2UI21|) are essential in order to investigate the thermodynamic limit properties 



for a fixed density N/L. They have the merit of being derived under mild assumptions for the 
two-body collision amplitudes and therefore with an enormous potential of applicability. 

In conclusion, we have solved exactly a system of spin-s particles that interact via arbitrary 
£7(1) factorizable ^-matrix. The structure of both eigenvectors and the eigenvalues was derived 
solely from the Yang-Baxter relations ()2I9|) and the unitarity property We believe that our 
formula (|13|) is capable to accommodate the solution of other integrable models possessing extra 
U(l) symmetries other than that already discussed. In these cases, previous experience [22] 



suggests that the recursive structure of (|T3j) will be preserved but now with i/)k(\i, \j 2 , . . . , Xj k ) 
behaving as a vector-function while B(A,/i) as an underlying factorizable S'-matrix. If this 
proposal turns out to be feasible in the future one would be able to solve exactly several 
distinct families of integrable models from a rather unified point of view. 
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